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A number of diverse bulk properties of the zincblende and wurtzite III-V nitrides A1N, GaN, 
and InN, are predicted from first principles within density functional theory using the plane-wave 
ultrasoft pseudopotential method, within both the LDA (local density) and GGA (generalized gra- 
dient) approximations to the exchange-correlation functional. Besides structure and cohesion, we 
study formation enthalpies (a key ingredient in predicting defect solubilities and surface stability), 
spontaneous polarizations and piezoelectric constants (central parameters for nanostructure model- 
ing), and elastic constants. Our study bears out the relative merits of the two density functional 
approaches in describing diverse properties of the III-V nitrides (and of the parent species N2, Al, 
Ga, and In), and leads us to conclude that the GGA approximation, associated with high-accuracy 
techniques such as multiprojector ultrasoft pseudopotentials or modern all-electron methods, is to 
be preferred in the study of III-V nitrides. 



I. INTRODUCTION AND METHOD 

The III-V nitride semiconductors A1N, GaN, and InN, 
and their alloys-are by now well established as a strategic 
material systema for applications in high-frequency opto- 
electronics (LED and LASER) , and high-power electron- 
ics (e.g. HEMTs). Most of their potential in these fields 
is due, respectively, to the large tunability of band gaps 
with alloy composition (in principle, 1.9 to 6.2 eV), and 
to their high peak and saturation drift velocity, coupled 
with polarization-induced effects allowing for the realiza- 
tion of high-density low-dimensional charge gases.H 

Nitrides physics posed a number of puzzles to (and 
profited considerably from) ab initio studies of various 
properties and subsystems, ranging from SAirfacesa to 
defectsjj and polarization-related properties.El Heralding 
the unusual nature of these materials, the standard study 
of the structural properties of bulk materials was a source 
of surprises in early-studies. For instance, since some of 
the earliest papersjj quite unusually for III-V semicon- 
ductors, the semicore 3c? electrons of Ga were found to 
behave as valence electrons and to be essential to describe 
accurately structural properties. 

A major source of uncertainty, both technical and ideo- 
logical in nature, in density-functional theory (DFT) cal- 
culations is the choice of the exchange-correlation func- 
tional. While LDA is used most commonly, the gen- 
eralized gradient approximation, or GGA, has become 
a close competitor in recent years. In this work, we 
study the effects of using either LDA or GGA in the 
prediction of the properties of III-V nitrides. Similar, 
calculations have been performed only once previously,Ll 
and were concerned with structural and cohesive prop- 
erties. In this paper, we add several new aspects to 



this theme. First, we iise ultrasoft pseudopotentials, 
which should in principleo improve over norm-conserving 
potentials .0 Second, we carefully calculate formation en- 
thalpies, which are the cornerstone of predictions on non- 
stoichiometric systems relevant to surface reconstruction 
and impurity solubility. This calculation requires us to 
study the metallic phases of Al, Ga, and In, the N2 
molecule, and solid nitrogen (a molecular solid compris- 
ing N2 dimers on an hep lattice). Third, we evaluate the 
spontaneous polarization and the piezoelectric constants 
of the wurtzite phasea in both the GGA and LDA. We 
find that these quantities are moderately affected by the 
choice of exchange-correlation, unlike most others prop- 
erties. Fourth, we evaluate a subset of the elastic con- 
stants in LDA and GGA. 

The calculations have been done using VASP (Vienna 
Ab-initio Simulation Package) ,E£l which implements the 
DFT scheme within both the LDA and GGA approx- 
imations: we adopted the well established PWfll ver- 
sion of the GGAU and Ceperley- Alder LOAM Ultra- 
soft pseudopotentials^ describe the electron-ion interac- 
tion. As usual, the potentials provided with VASP are 
generated for the free atom using the appropriate (LDA 
or GGA) functional. The pseudopotentials for Ga and 
In include respectively the semicore 3d and Ad states in 
the valence. A plane wave basis is used to expand the 
wave- functions. We use a cut off of 350 eV, which is suf- 
ficient to fully converge all properties of relevance. For 
k-space summation, we use at least a Monkhorst-Pack 
(888) grid, except for the N2 molecule. Lattice constants 
and internal parameters are calculated using standard to- 
tal energy calculations. Polarizations are obtained using 
the Berry phase approachtS as in previous work.B Cohe- 
sive energies are calculated relative to spin-polarized free 
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atoms. The formation enthalpies Ai/xN per atom pair 
of the XN crystals are calculated as 
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where Exn is the total energy per atom pair of the com- 
pound XN, Ex the energy per atom of bulk X = Al, Ga, 
and In, and E^ is the energy per N atom in the N2 dimer 
or the condensed N2 phase. 



II. PARENT SPECIES 
A. Nitrogen: molecule and solid 

The nitrogen dimer is studied in artificial periodic con- 
ditions in a cubic box of side 10 A, using the T point for 
k-summation. The results, listed in Table |, agree well 
with other LDA and GGA calculations. GGA shows an 
overall better agreement with experiment. The binding 
energy is evaluated including the spin-polarization energy 
of the N atom (-2.89 eV), calculated with a locaLspin 
density all-electron scalar-relativistic atomic codeJiil 

Solid nitrogen is a condensate of N2 molecules. We 
consider the stable phase, with vertically-oriented N2 
molecules centered at the lattice points of a close packed 
hexagonal lattice. The (888) grid is used for k-space sum- 
mation. We compare our results with experimental data 
from Ref. [l5] in Table || From the data in Tables | 
and [n|, the binding energy per molecule in the condensed 
phase is 0.328 eV in the LDA and 0.143 eV in the GGA. 
While performing well as to the in-plane lattice constant 
(about ±0.5% relative deviation), both functionals fail to 
some extent with the axial lattice parameter: LDA un- 
derestimates it strongly (~20%) and GGA overestimates 
it (~9%). The vertical center-to-center intermolecular 
distances are 3.33 A experimentally, 3.64 A in GGA, 
and 2.68 A in LDA. This system is indeed a severe test 
for both functionals because of its weak dipolar binding. 
GGA performs slightly better, as expected. The binding 
of the N2 system, already extremely large in reality, is 
overestimated appreciably in both approaches. 



B. Bulk Al, Ga, and In 

Metallic Al, Ga and In are a necessary ingredient to 
calculate formation enthalpies. Al is a good free electron 
metal, and has the fee structure. Ga is a mixed-bonding 
marginal metal (see e.g. Ref. |l6|). At ambient conditions, 
its stable phase is a dimerized structure known as a-Ga, 
a face-centered orthorombic lattice with crystallographic 
vectors ai = ax, a 2 — ^(by + cz), a 3 = |(— by + ci) and 
eight atoms per primitive cell, whose positions ar e d e- . 
fined by two additional internal parameters u and t>.li3'L£l 
Indium crystallizes in the monoatomic body-centered 
tetragonal lattice, with lattice constants a and c. and 
primitive vectors ai = ax, a 2 = ay, a3 = ^(ax+ay+cz). 



We report our results for Al in Table III , for Ga in 
Table |v|, and for In in Table [vj In these calculations we 
use a cutoff of 350 eV, an (888) k-space mesh for Al and 
In, and a (12,12,12) mesh for Ga. The cohesive energy 
includes the atomic spin-polarization energy (-0.136 eV 
for Al, -0.134 eV for Ga, and -0.117 eV for In). As usual 
GGA improves somewhat the lattice constant and bind- 
ing energy. For Ga, both approaches are off by about 
the same amount in opposite directions for a. In both 
cases, axial ratios and internal parameters are excellent. 
Our LDA results improve somewhat over those ofRef.|l|, 
presumably because of the explicit treatment of 3d elec- 
trons. For In, LDA and GGA are again off the mark by 
equal and opposite amounts for a. The LDA axial ratio 
is slightly better than GGA's. In short, the usual trend 
is obtained of expanded and softer lattice as produced by 
GGA compared to LDA. If one is forced to choose, GGA 
generally performs better, especially in terms of cohesive 
energies. In any case, the deviations typically are below 
± 1%, so both approaches are quite legitimate. 



III. THE NITRIDES 

Binary III-V nitrides occur in nature in the wurtzite 
structure (the (3 phase). Zinc-blende nitrides (the a 
phase) have a slightly higher energy. It is possible to 
grow epitaxially, e.g., a— GaN on cubic substrates. We 



first analyze zinc-blende (Sec. Ill A), then wurtzite (Sec. 
Ill B|) . Our results are compared with those of Ref. [7], 



where numerous other theoretical values are provided. 



A. Zinc-blende A1N, GaN, InN 

For zinc-blende nitrides we used the usual 350 eV cutoff 
and (888) k-grid. To estimate the cohesive energy, we use 
the atomic spin-polarizations ind icat ed pr eviou sly. Our 
results are reported in Tables VI, VII, and VIII, for A1N, 
GaN, and InN respectively. The results confirm the by 



now usual behavior of GGA vs LDA, consisting in a soft- 
ening of the lattice, which improves lattice constant and 
binding energy, but worsens slightly the bulk modulus. 
Comparing the cohesive energy with that of the wurtzite 
phase as discussed below, we find that zinc-blende is dis- 
favored over wurtzite. 



B. Wurtzite A1N, GaN, InN 

Wurtzite is a hexagonal close-packed lattice, compris- 
ing vertically-oriented X-N units at the lattice sites. The 
basal lattice parameter is a, the axial lattice parameter is 
c. The interatomic distance in the basic unit is described 
by an internal parameter u expressed in units of the axial 
ratio c/a. The ideal (i.e. for packed hard spheres) values 
of the axial ratio and internal parameter are respectively 
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c/a=y/8/3 and u — 3/8. The crystallographic vectors of 
wurtzite are a=a (1/2, VH/2, 0), b=a (1/2, -a/3/2, 0), 
and c=a (0, 0, c/a). The cartesian coordinates of the 
basis atoms are (0, 0, 0), (0, 0, uc), a (1/2, v3/6, c/2o), 
and a (1/2, a/3/6, [tt+1/2] c/a). 

Our results are reported in Tables [lX|, [x|, and XI , for 
A1N, GaN, and InN respectively. For comparison, exper- 
imental data, and the results of Ref. |?| are also listed. As 
to structure, in all cases both the axial ratio and the in- 
ternal parameter are non-ideal. Deviation from ideality 
increases from GaN to InN to A1N. As usual GGA im- 
proves considerably lattice constant and binding energy, 
at the cost of a slight overestimate of the axial ratio. 
The internal parameter u (alias the axial bond length) is 
well reproduced in all the various combination of materi- 
als and approximations. The present GGA calculations 
produce lattice constants and internal parameters with 
maximum deviations from experiment below 0.3 % for 
A1N, 0.9 % for GaN, and 1.7 % for InN (+1.3% for a 
and +1.7% for c). In this respect, these are probably the 
best DFT pseudopotential results so far for these ma- 
terials. The improvement over previous GGA results is 
to be attributed pto the use of ultrasoft, multiprojector 
pseudopotentialsu By the same token, it is quite likely 
that all-electron calculations using the same GGA pa- 
rameterization may improve the agreement further, es- 
pecially for InN. 

Calculated cohesive energies overestimate, as usual, 
the experimental value. GGA corrects in part the LDA 
overbinding, and exhibits better agreement. Compar- 
ing the cohesive energies of the zinc-blende and wurtzite 
phases, as already mentioned, we find wurtzite to be en- 
ergetically favored over zincblende. The predicted differ- 
ence per atom pair between the two phases is 189 meV 
(LDA) and 164 meV (GGA) for A1N, 17 meV (LDA) and 
16 meV (GGA) for GaN, 17 meV (LDA) and 15 meV 
(GGA) for InN. 

Good results are also achieved for the formation en- 
thalpies. These were obtained using the energy per N 
atom in the solid-N2 phase: if the N2 molecule is as- 
sumed as reference instead, one half of the binding en- 
ergy of solid-N 2 (i.e. 0.164 eV in LDA, and 0.071 eV in 
GGA) must be added to the values in the Tables. The 
reported values are referred to an atom pair. The avail- 
able formation enthalpy measurements for the nitrides 
are not very recent Jl3-t3 and their spread is considerable 
(this is presumably also the case for cohesive energies, 
which are measured as heats of vaporization) . The cal- 
culated GGA formation enthalpies are in general agree- 
ment with experiment for A1N and GaN. For InN, GGA 
overcorrects the LDA overbinding and gives a positive 
value. This problem is due to InN itself, and not to In or 
N parent phases: indeed, quite unusually, the calculated 
cohesive energy underestimates the experimental value. 
We are not aware of other formation enthalpy calcula- 
tions for InN. The issue is open to further investigation, 
especially by all-electron methods. 



In Table XII, we report for each of the nitrides the 



spontaneous polarization in the equilibrium structure, 
the dynamical effective charges, the piezoelectric con- 
stants, and a subset of elastic constants relevant to 
symmetry-conserving strains. The reason for collecting 
these data in one Table is that they provide an almost 
self-contained set of input data for the simulation of 
nanostructures made of wurtzite nitrides. The only ad- 
ditional data needed are the static dielectric constants, 
which were reported elsewhere.tj In the last column we 
report the,pjpper piezoelectric constant e^. As discussed 
recentlyJ23'EZI this value should be compared with experi- 
ments involving current flow across the sample, whereas 
the "improper" constant e%i is relevant to systems in de- 
polarizing fields such as nitride nanostructures .LI 

It is not infrequent to hear the incorrect statement that 
the spontaneous polarization is non vanishing in wurtzite 
because of structural non- ideality. In actual fact, a noiu 
vanishing polarization is allowed on symmetry groundsel 
in the ideal wurtzite structure as well. Indeed, we find 
that the calculated Berry-phase polarization in the ideal 
structure is indeed non-zero: -0.032 C/m 2 in A1N, -0.018 
C/m 2 in GaN, and -0.017 C/m 2 in InN. These values are 
smaller (by a factor of 2 t o 3) than the actual ones for 
non-ideal structures (Table XII). This confirms the intu- 
itive idea that non-ideality, and especially changes in u, 
can increase polarization substantially, and indicates that 
an accurate determination of the structure is mandatory 
to obtain reliable polarization values EO 

Theoretical predictions on polarization properties were 
shown to compare quite favorably with experimental ev- 
idences in various papers (see e.g. Refs. || |, and ||). It 
should be noted, however, that the link between polariza- 
tion and the observed quantities, typically optical shifts 
or densities of mobile charge, is rather indirect and af- 
fected by uncertainties due to issues of nanostructure de- 
sign, material quality, and reverse modeling. Thus, com- 
parison with experiment does not yet allow a clear-cut 
evaluation of the. performance of LDA vs GGA. The re- 
cently reportedcJ non-linear behavior of the polarization 
in nitride alloys is an additional source of uncertainty. 

The LDA elastic constants are in fair agreement with 
those of Wright.cil The GGA constants are smaller, as 
is to be expected given the general tendency of GGA to 
produce a softer lattice. According to elasticity theory, 
the axial strain induced in wurtzite by an in-plane (e.g. 
epitaxial) strain t\ is 63 = — 2 e\ C^i / = Re%. The 
quantity R is thus relevant to epitaxial nitride systems, 
and it is reported in Table |XII| . Several experimental 
data for the elastic constants and R are compiled in Ref. 
p4| . The considerable spread of those data does not al- 
low a definite conclusion about whether GGA produces a 
systematically improved agreement with experiment over 
LDA in this respect. 
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IV. SUMMARY 

In conclusion, GGA and LDA calculations for III-V 
nitrides suggest an overall improvement of the predicted 
properties in the former approximation. In particular, 
the structural parameters are extremely accurate, with 
deviations from experiment below 0.3 % for A1N, 0.9 % 
for GaN, and 1.7 % for InN. Cohesive energies and for- 
mation enthalpies are in fair to excellent agreement with 
experiment; the only clear-cut failure (in the GGA ap- 
proximation) is the positive formation enthalpy of InN. 
Elastic properties follow the expected trends of GGA vs 
LDA behavior; due to uncertainties in the experimental 
data, comparison with experiment does not provide defi- 
nite support to one or the other approximation. Polariza- 
tion properties are moderately sensitive to the exchange- 
correlation functional, as long as the latter predicts the 
correct structure (especially, the correct internal param- 
eter u). For these properties too, comparison with exper- 
iment is indirect and affected by many sources of uncer- 
tainty, and does not support one or the other approach. 
Concerning the cohesion and structure of parent species 
(N 2 , Ga, Al, In), only in the case of condensed N 2 do we 
find major discrepancies with experiment. In the light 
of the present results, our conclusion is that the GGA 
approximation should be preferred in density-functional 
studies of III-V nitrides. Of course, the choice of one 
or the other approximation may depend on the specific 
problem being addressed. 
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TABLE I. Bond length, vibrational frequency, and binding 
energy of the N2 dimer. 





d(A) 


lo (THz) 


E„ (cV) 


LDA 


1.107 


464.3 


-11.332 


GGA 


1.113 


442.8 


-10.558 


Experiment 51 


1.10 


444.8 


-9.9 


a Ref. | 








TABLE II. Structural paramet' 


ers and binding > 


energy per 


molecule of hexagonal solid N2 








a (A) 


c/a 


E„ (eV) 


LDA (present) 


4.0205 


1.3311 


-11.660 


GGA (present) 


4.0633 


1.7929 


-10.701 


Experiment 3 


4.039 


1.6514 




a Ref. 








TABLE III. 


Lattice constant, 


binding energy, 


and bulk 


modulus of bulk fee Al. 








a (A) 


E b (eV) 


B(MBar) 


LDA (present) 


3.9809 


-4.064 


0.766 


GGA (present) 


4.0491 


-3.561 


0.689 


Experiment 3 


4.05 


-3.39 


0.773 



a Ref. 



TABLE IV. Lattice constant, binding energy, axial ratios, 
and internal parameter (units of c) of a-Ga. 





a (A) 


b/a 


c/a 


u 


V 


E„ (eV) 


LDA (present) 


4.4365 


0.9985 


1.6856 


0.0816 


0.1577 


-3.484 


LDA a 


4.377 


0.994 


1.688 


0.0803 


0.1567 




GGA (present) 


4.5962 


0.9917 


1.6961 


0.0834 


0.1559 


-2.796 


Experiment 3 


4.51 


1.0013 


1.695 


0.0785 


0.1525 





a Ref. 



TABLE VI. 


Lattice constant, 


binding energ 


;y, and forma- 


tion enthalpy AH of zinc-blende A1N. 






a (A) 


E 6 (eV) 


A ZJ /~"\A 

£\Jrt ^evj 


LDA (present) 


4.332 


-13.347 


-3.449 


LDA a 


4.310 


-13.242 




GGA (present) 


4.390 


-11.907 


-2.975 


GGA a 


4.394 


-11.361 




Experiment 1 " 


4.37 






a Ref. U 








b Ref. ]1§ 








TABLE VII. 


Lattice constant, binding energy, and forma- 


tion enthalpy AH of zinc-blende GaN. 






a (A) 


E 6 (eV) 


AH (eV) 


LDA (present) 


4.446 


-10.982 


— i.Doy 


LDA a 


4.518 


-10.179 




LDA b 


4.466 


-10.880 




GGA (present) 


4.538 


-9.249 


-1.102 


GGA a 


4.590 


-8.253 




Experiment 


4.519 






a Ref. | 








b Ref. |L 








c Ref. y 








TABLE VIII 


. Lattice constant, binding ener 


gy, and forma- 


tion enthalpy AH of zinc-blende InN. 






a (A) 


Ef, (eV) 


AH (eV) 


LDA (present) 


4.964 


-9.232 


-0.282 


LDA a 


5.004 


-8.676 




GGA (present) 


5.067 


-7.680 


0.140 


GGA a 


5.109 


-6.855 




Experiment 13 


4.98 







a Ref. M 
b Ref. jij] 



TABLE V. Lattice constant, axial ratio, and binding en- 
ergy of bulk In. 





a (A) 


c/a 


E b (eV) 


LDA (present) 


3.1861 


1.5348 


-3.116 


GGA (present) 


3.2958 


1.5448 


-2.470 


Experiment 


3.244 a 


1.5222 a 





'•Ref. [l| 



TABLE IX. Lattice constant, axial ratio, internal parame- 
ter, and formation enthalpy of wurtzite A1N. 



a (A) c/a u 



E 6 (eV) AH (eV) 



LDA (present) 3.0798 1.5995 0.3821 -13.536 -3.642 

LDA a 3.057 1.617 0.3802 -13.286 

GGA (present) 3.1095 1.6060 0.3819 -12.071 -3.142 

GGA a 3.113 1.6193 0.3798 -11.403 

Experiment 11 3.1106 1.6008 0.3821 c -11.669 a -3.13 d -2.53 e 
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TABLE X. Lattice constant, axial ratio, internal parame- 
ter, and formation enthalpy of wurtzite GaN. 



a (A) c/a 



E(ev) AH (eV) 



LDA (present) 3.1461 1.6273 0.3768 -10.999 -1.685 

LDA a 3.193 1.634 0.376 -10.187 

GGA (present) 3.1986 1.6339 0.3772 -9.265 -1.118 

GGA a 3.245 1.632 0.3762 -8.265 

Experiment 13 3.1890 1.6263 0.377 -9.058 a -1.08 c -1.91 d 




TABLE XL Lattice constant, axial ratio, internal parame- 
ter, and formation enthalpy of wurtzite InN. 



(A) c/a 



E b (eV) AH (eV) 



LDA (present) 3.5218 1.6121 0.3791 

LDA a 3.544 1.626 0.377 

GGA (present) 3.5848 1.6180 0.37929 

GGA a 3.614 1.628 0.377 

Experiment 13 3.538 1.6119 

a Ref. @_ 
b Ref. 
c Ref. 
d Ref. 
°Ref. 



-9.249 
-8.694 
-7.695 
-6.872 
-7.970 c 



-0.303 



0.125 



-0.21 d -1.36 s 



TABLE XII. Spontaneous polarization (C/m 2 ), piezoelec- 
tric constants (C/m 2 ), dynamical charges, elastic constants 
(GPa), and the ratio R=— 2 C31/C33 (see text) of wurtzite ni- 
trides, as obtained in the LDA and GGA approximation. The 
last column reports the proper esi piezoelectric constant. 





P 


Z* 


e33 


e3i 


C33 


C31 


R 


e 31 


A1N 


















LDA 


-0.100 


2.652 


1.80 


-0.64 


384 


111 


-0.578 


-0.74 


LDA a 










373 


108 


-0.579 




GGA 


-0.090 


2.653 


1.50 


-0.53 


377 


94 


-0.499 


-0.62 


GaN 


















LDA 


-0.032 


2.51 


0.85 


-0.44 


415 


83 


-0.400 


-0.47 


LDA a 










405 


103 


-0.508 




GGA 


-0.034 


2.67 


0.66 


-0.34 


354 


68 


-0.384 


-0.37 


InN 


















LDA 


-0.041 


3.045 


1.09 


-0.52 


233 


88 


-0.755 


-0.56 


LDA a 










224 


92 


-0.821 




GGA 


-0.042 


3.105 


0.81 


-0.41 


205 


70 


-0.683 


-0.45 



a Ref. 24 
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